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A novel description of Josephson vortices (JVs) crossed by the pancake vortices (PVs) is proposed 
on the basis of the anisotropic London equations. The field distribution of a JV and its energy 
have been calculated for both dense (a < Aj) and dilute (a > A./) PV lattices with distance a 
between PVs, and the nonlinear JV core size Xj. It is shown that the "shifted" PV lattice (PVs 
displaced mainly along JVs in the crossing vortex lattice structure), formed in high out-of- plane 
magnetic fields B z > $o/7 2 s 2 [A.E. Koshelev, Phys. Rev. Lett. 83, 187 (1999)], transforms into 
the PV lattice "trapped" by the JV sublattice at a certain field, lower than <E>o/7 2 s 2 , where $o is 
the flux quantum, 7 is the anisotropy parameter and s is the distance between CUO2 planes. With 
further decreasing B z , the free energy of the crossing vortex lattice structure ( PV and JV sublattices 
coexist separately) can exceed the free energy of the tilted lattice (common PV-JV vortex structure) 
in the case of 7s < \ a t with the in-plane penetration depth A a f, if the low (B x < 7 ( l > o/A 2 i) ) or high 
{B x > <&o/7s 2 ) in-plane magnetic field is applied. It means that the crossing vortex structure 
is realized in the intermediate field orientations, while the tilted vortex lattice can exist if the 
magnetic field is aligned near the c-axis and the a&-plane as well. In the intermediate in-plane 
fields 7$ /A 2 b < B x < $ /ys 2 , the crossing vortex structure with the "trapped" PV sublattice 
seems to settle in until the lock-in transition occurs since this structure has the lower energy with 
respect to the tilted vortex structure in the magnetic field H oriented near the afo-plane. The 
recent experimental results concerning the vortex lattice melting transition and transitions in the 
vortex solid phase in Bi2Sr2CaCu2 08+i single crystals are discussed in the context of the presented 
theoretical model. 
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The mixed state of high temperature supercondjucr 
tors is complex and rich with various vortex phasesErH. 
Besides the vortex lattice described by 3D anisotropic 
Ginzburg-Landau model, the new types of vortex struc- 
tures can occur within the large part of the phase diagram 
of the mixed state where the coherence length along the 
c-axis is smaller than the distance between CuC>2 planes. 
In such a case, the magnetic field, aligned with the c- 
axis, penetrates a superconductor in the form of quasi 
two-dimensional pancake vortices (PVs)B while the field 
applied parallel to the a6-plane generates Josenhson vor- 
tices (JVs) in the layers between Cu02 planeaTEl In mag- 
netic fields tilted with respect to the c-axis, PVs and JVs 
can form a common tilted latticeu or exist separately as 
a crossing (combined) latticeoQ. The tilted lattice rep- 
resents the inclined PVs stacks in fields applied close to 
the c-axis while, at higher angles, the pieces of JVs link- 
ing PVs are developedtffi. The crossing lattice is another 
structure containing both a PV stack sublattice and a JV 
sublattice which coexist separately. 

The vortex-solid phase diagram in the tilted mag- 
netic fields was first proposed by Bulaevskii, Ledvij, and 
KogarJj. According to their model, which does not take 
into account the interaction between PV and JV sub- 
lattices in the crossing lattice structure, the tilted lat- 
tice is formed for all orientations of the magnetic field 
until the lock-in transitionEl occurs if the in-plane Lon- 
don penetration depth A a & is larger than the Josephson 
vortex core with size 7s (7 is the anisotropy parame- 



ter and s is the distance between CuC>2 planes). In 
the opposite limit, 7s > A a b, the tilted lattice trans- 
forms into the crossing lattice (as the magnetic field is 
inclined away from the c-axis) _at a certain angle before 
the lock-in transition happensO. Later, the possibility 
of the coexistence of two vortex sublattices with differ- 
ent orientations was analyzed numerically by comparing 
the free energy of such system with the free energy of 
mono-oriented tilted vortex lattice at different field ori- 
entations and different absolute values of the external 
magnetic field for the case of 3D anisotropic (London 
model) superco^uctorsoilil as wellras layered (Lawrence- 
Doniach-jaaadeJiia) superconductorstS. According to that 
analysiaHllij performed for 7 = 50 — 160, the crossing 
lattice can be energetically preferable in the quite low 
magnetic fields (B = y 1 B 2 Z + B 2 < $o/\ 2 ab ) in the in- 
termediate field orientations < 6% < 8 < 62 < tt/2 
with 9 = arctan(i? x /_B z ) (B z and B x are the field compo- 
nent along the c-axis and parallel to the afc-plane, respec- 
tively). However, the interaction of two coexistiaej-Mortex 
sublattices was apt considered in those worksmliLj. Re- 
cently, KoshelevQ has studied the case of extremely ani- 
sotropic superconductors 7s \ a b and has shown that 
the crossing lattice can occupy substantially larger region 
of the vortex lattice phase diagram in the oblique fields 
due to the renormalization of the JV energy £j through 
the interaction of a Josephson vortex and the PV sublat- 
tice. In additian_such interaction leads to the attraction 
of PVs to JVaL-3'Ll at low out-of-plane magnetic fields B z 
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(the some sort of pinning effect). This pinning may in- 
duce transitions between different substructures of the 
crossing lattice structure. However, there is still no the- 
oretical investigation how PV sublattice can influence on 
the JV lattice in the crossing vortex structure in the case 
of moderate anisotropic superconductors with 7s < X a b- 
In this regard, the phase diagrams of the vortex lattice 
for strongly anisotropic layered superconductors should 
be reconsidered (at least for the case of 7s < A a b) by 
taking into account the renormalization of £j and the 
pinning of PVs by J Vs. 

The vortex structures in highly anisotropic layered su- 
perconductors are usually studied on the basis of the 
nonlinear discrete Lawrence-Doniach modelll3, but this 
model is quite complex and the detailed analysis of the 
vortex system is complicated. On the other hand, the 
layerness of superconductors can be ignored on scales 
larger than the size of the nonlinear Josephson vortex 
core. Therefore, the linear anisotropic London model 
could be applied for a study of the vortex crossing lattice 
outside JV cores. In this paper we introduce the extended 
London theory which allows to describe the crossing lat- 
tice as well as to calculate the energy £j and the field 
distribution of JV in the presence of the crossed PV sub- 
lattice. It is shown that with decreasing the perpendicu- 
lar magnetic field the pancake sublattice transforms from 
the "shifted" sublattice characterized by one component 
PV displacement along JVs to the "trapped" sublattice 
where JVs are occupied by PV rows. The comparison of 
the free energies of the tilted lattice and the crossing vor- 
tex structure for the case X a b > 7s indicates that in low 
(Bx < l^o/Kb) and high (B x > <$> /js 2 ) in-plane fields, 
the tilted lattice can exist if the vector B = B x e x +B z e z 
is directed close to the c-axis as well as near the afr-plane, 
while the crossing lattice is realized in the fields oriented 
far enough from the crystal symmetry axes. Further- 
more, in the intermediate in-plane fields the tilted vortex 
lattice exists only at the magnetic field orientations near 
the c-axis whereas crossing vortex structure settles in the 
wide angular range until the lock-in transition happens. 

This paper is organized as follows. The general equa- 
tions for the magnetic field distribution and the energy 
of Josephson vortex in the presence of the pancake lat- 



The dense pancake lat- 
It is shown that, in the 



tice are derived in section 
tice is studied in section ft 
limit 7s > A q6 and $0/(75? < B z < ^o7 2 s 2 /X ab , our 
model reproduces the results which were earlier obtained 
by KoshelevQ, while the shear deformation of the PV lat- 
tice significantly renormalizcs the JV energy at the higher 
out-of-plane fields. Section [II is devoted to the dilute PV 
lattice. It is described how the novel vortex substructure 
with the PV lattice "trapped" by the JV lattice can be 
realized at low B z . The phase diagram of the vortex-solid 
phase in the tilted magnetic fields is considered for the 
case of X a b > 7s in section IV while the recent experi- 
mental results are discussed in section [v|. 



I. JOSEPHSON VORTEX IN THE PRESENCE OF 
PANCAKE VORTEX LATTICE: GENERAL 
EQUATIONS 

We consider a Josephson vortex crossed with the pan- 
cake lattice in the framework of the modified London 
model. On scales which are much larger than both the 
distance between Cu02 planes and the in-plane coher- 
ence length £ a {,, the pancake vortex stack could be con- 
sidered as an ordinary vortex line at temperatuces-pgnif- 
icantly lower than the evaporation temperatureuHj. The 
same approach can be also used for the description of 
the Josephson vortex far from the nonlinear core. The 
JV current acts on PVs through the Lorentz force causing 
their displacements along JV, which can be interpreted as 
a local inclination of the PV lines away from the c-axis. 
In turn, the local tilt of the PV stacks induces an ad- 
ditional current along the c-axis which redistributes the 
"bare" JV field. Such physical picture can be described 
with one-component PV displacement u = (u, 0, 0) which 
does not depend on the x-coordinate (Fig. 1). The free 
energy functional Tpj can be written as 
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Tpj=—\ d J R(h; + Vxh p AVxhp 
+ h,/ 2 + VxhjAVxhj + 2h p hj + 2Vxh„AVxhj ) , 



(1) 

where h p and hj are the magnetic fields of PV lines 

and JV, respectively, and A is the penetration-depth ten- 
sor, V = (d/dx,d/dy,d/dz). In the considered coordi- 
nate system the tensor has only the diagonal components 



A 



A„ 



A 



b , A zz — X c with anisotropic penetration 
depths X a b and A c . The field h p is determined by the dis- 
placement u„of PVs through the London equation (see, 
for instance,B) 

hp+\7x ^AVxhp^j = § ^Jdz 

du(Yi,z)e x 



di 



6(r-K z {S)-u(Y l ,S)e x ), (2) 



where Rj(z) = (Xi,Yi,z) is the equilibrium position of 
the z-th PV line, r = (x,y,z) while e z and e x are unit 
vectors along the z and x axes, respectively. (Here we 
have accepted that the parametric equation 17 (z) (with 
parameter z) describing the i-th vortex line takes the 
form of x{z) = Xi + Ui(z), y(z) = Yi, z(z) = z.) The 
in-plane coordinates of the unshifted lines Xi and Yi are 
expressed through the distances a and b (see Fig. I b) 
between PVs and PV rows as Xi = al/2 + aj and Yi = bl 
with integer I and j. In our approach, the field of the 
Josephson vortex also obeys the London equation 



A 



d 2 h r 

' dz 2 



A; 



d 2 h.j 
dy 2 



*o8(y)6(z), 



(3) 
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FIG. 1. The JV crossed by PV stacks which are shifted by 
the currents induced by JV: a) 3D sketch depicts the defor- 
mation of the PV lattice in the different CuC>2 planes, b) 2D 
sketches show the deformation of the PV lattice in a CuCh 
plane for the dense (left) and dilute (right) PV lattices. The 
filled circles correspond to the unshifted PVs while the open 
ones represent the PVs shifted due to the interaction with JV 
currents. The shaded area images the nonlinear JV core re- 
gion. The dashed-dotted lines mark the rows of the unshifted 
PVs while the dashed curve shows the deviation of PVs from 
these rows. 



where <5-functions should be smoothed on a scale of the 
Josephson vortex core. The JV core size along the z-axis 
is fixed by the interlayer distance s, while the core length 
along the y-direction is limited by the condition that the 
current along the c-axis can not exceed the maximum 
interlayer current j c ~ c$o/(8ir 2 A 2 s). In the presence 
of PVs, the current across the layers consists of both 
the current of JV itself and the current born by the local 
inclination of PV lines. Therefore, the core size along the 
y-axis A j can be renormalized in the presence of PVs and 
should be calculated self-consistently (see next section). 
Furthermore, the space variable y could be replaced by 
y — yo in the argument of the (^-function with < yo < b 
since the PV lattice can be arbitrary shifted from the 
center of JV. However, we take yo — which corresponds 
to the energetically more preferable positionQ. 

Next, in order to find the distribution of the magnetic 
field in the vortex system and the energy of JV, we will 
minimize the free energy (^) as a functional of the dis- 
placement u. The fields h p and hj can be obtained using 
equations (00) with the displacement u fixed by the min- 
imization of (gj). Then, the energy of JV, £j, defined as 
the difference of the free energies (|l|) with and without 
JV, will be derived. This energy includes the self energy 
of JV and the change of the free energy of the PV lattice 
born by the interaction with JV. We will use the elastic 
approximation, i.e., the free energy ([!]) and the magnetic 
field of PVs (0) will be expanded up to the second order 



Using the integral representation of 5-function, the 
equation can be rewritten as 



hp+Vx ^AVxhp^ =®o^2f d 

x f _^!3_ e iqr -iqRi(i) ( M 

J (2^ 6 \ z( 



iq x Ui(z) - -3 x Ui(2)) 



duAz). 
az 



(4) 



The field h p of PV lines changes on different space scales. 
The first scale is determined by the characteristic gradi- 
ent of the displacement u(y, z) and usually is much larger 
than the distance a between PVs. The second scale is de- 
fined by the discreteness of the PV lattice and it is about 
a. To separate the contribution to the free energy from 
these scales, we introduce the Fourier variables u(k y , k z ): 



u(Yi,z) = 



J -rib 2n J 2n 



and 



i(k y ,k z ) = b}^ / dzu(Yi, 



z)e 



—i(k y Yi-\-k z z) 



(5) 



(6) 



where the domain of variation of k z is restricted 
by the inequality \k z \ < 1/s born by the layerness 
of the system. Substituting (g) into (J|) and using 
the well-known equality dz exp mk — q)Ri(z)) = 

(2tt) 3 (B z /$ ) J2 q $(k - q - Q), where Q = (Q x ,Q y , 0) 
are the vectors of the reciprocal lattice (Q x = 
27T77i/a, Q y — 7r(2n + m)/b with integer m and n), one 
gets the expansion of the field of PVs in series with re- 
spect to the displacement u: 



h p = h<°> + h« [u] + h p 2 ) [u] , rip = n<°> 
hW + VxAVxh^ = nW 



P 

where 

v, (0) 



(7) 



i 

Q J 

x (e z iQ x + e x ik z )u{k y ,k z )e- lQ * x e^ k y- Q ^e lk *\ 

n ( 2 ) = B z Y [ 
p ^ J 16?r 4 



dkydk z dk y dk z 

-u(k)u(k ) 



x I ~^Q x ^z Q x k z e 3 



e i(k z +k' z )z e -iQxX e i(k v +k' v -Q v )y 



(8) 
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The term of h p with Q = corresponds to the continuous 
approximation and varies on the large scale, while terms 
with Q^O are related to the field components changing 
on the scale of a. 

It is easy to see that only terms with Q x — will 
give contribution to the part of the free energy describ- 
ing the interaction between PVs and JV, because the 
field hj does not depend on x and all terms with Q x ^ 
vanish after integration over x. Therefore, it is con- 



venient to divide hp 1 -* and tl£' into two components 



(i) 



,(<9) 



h*, where h p 



(Q) 



and 



rip"" include summands with Q x ^ while h* and n* do 
not vary with x. Then, the free energy functional F cross , 
containing only terms dependent on the displacement u, 
can be introduced as F cross = Tpj — Tp — Tj, where 
Tp and Tj are the free energies of the unperturbed PV 
lattice and the "bare" JV, respectively. Using equations 
(0), we obtain the expression for F cross as 



Across = ^ / d 3 R (n««h(« + 2h(°)n p 2 ) 



1 

8^ 



d 3 R (h* p n* p + 2hjn* p ) 



(9) 



The first contribution comes from the terms with Q x ^ 
and depends only on the short-scale variations of h p . It 
is determined by the shear deformation and the tilt de- 
formation. The second part describes the interaction of 
PVs with JV and with the current generated by PVs 
along the y-axis. Using the equations (||), the free en- 
ergy F cross can be rewritten in term of Fourier variables 
u(ky, k z ): 



F 



1 f dk dk 

' ' " 5S = 2 / 471-2 ~ ( u ee(k y ) + Uu(ky,k z )) u(k)u(-k) 



47T 



E 



dkydk z , . 



Q*=0' 

f(k z , k y - Q y ) - (g z /2$ )ifc z u(-k) 
l + Xl b kl + Xl{k y -Q v f 

with the shear energy 



(10) 



66 




i + K b Ql + K b (ky-Qy) 2 ) 



i + K b (Ql + Ql) 



(ii) 



and the tilt energy 



Ql 



UM ^ Q 5o I 1 + W + X lbQl + K b (*v Qv) 2 

Ql 



l + K b Ql + K b {Qy-ky) 



1 + Xl b kl + XIQI + Xl(k y - 

2 1,2 



(K-K b )Qlk 



(1 + Xl b kl + Xl b Ql + Xl b (k y 



,) 2 ) 



(l + \lQl + \l{ky-QyY + Xl b kl) 



(12) 



The expressions for 1/44 and Uqq represent sums over the 
reciprocal lattice vectors with Q x ^ 0, while the sum- 
mation in the last term of (10) is performed only over 
the reciprocal lattice vectors with Q x = 0. The function 
/(q) in ( |l0| ) appears due to smoothing of S function in 
eq. (||) and can be evaluated as /(q) « 1 in the rectan- 
gular region \q z \ < 1/s, \q y \ < 1/Xj and / ~ outside 
that area. The summation in the expression ( |Ti"| ) for the 
shear elastic energy was done by BrandtEj in the limit 
k y <C ir/b: 

U e6 = C 66 k 2 y , (13) 

where the shear elastic modulus Cqq is expressed as 
C 66 = (B 2 $ )/(87rA a6 ) 2 for a = y/$ /B z < X ab , 



while C 66 = ^/irX ab /(6a )& /(4irX 2 ab ) 2 
a o > ^a6- The tilt energy was obtained i 



a /A ab ) for 



U * = ii$k r 1 1 



k z K bln 



K b 2 + K 2 Q 



£ ab 



(14) 



\l V^o 2 + (^ 2 /7 2 ) + A c - 2 , 
for k z > Kq — 2ir/b, while 



(47rA ah ) 4 



32tt 2 A 2 . 



2 

44 n z 



(15) 

for k z <C Kq. Performing summation over Q y in the 
second term of equation (|l^) (see Appendix A) , we finally 
obtain the free energy functional: 



dky dk 2, 
~2it^2tt 



1 



(U M + U 66 ) u(k)«(-k) 



+E± lkz y {kzj ky) ^( k ) _ ffc, J|- u (k) u (-k) ) I (i 6) 
where $ is defined by the equation 



y(k z ,k v ) 



2A c yr+Api 



cosh I + ^ 2 b k%b/X c J — cos k y b 



(17) 
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for k y < min(7r/6, 1/Aj) and k z < l/s while \& ~ out- 
side that rectangular area. In the case of small values of 
wave vector k (k y <C 7r/6 and k z <C 7/6), the discreteness 
of PV lattice is irrelevant and the function ^> coincides 
with the Fourrier image of the "bare" JV field, but \P is 
modified substantially for larger k y or k z . 

The minimization of the free energy functional ( |l6| ) 
determines the displacement u as 



scales larger than the distance between PV lines even 
near the JV core. Thus, the continuous approximation is 
applicable in the whole space. In this case \k y \ < 1/Xj < 
n/b and, therefore, the cosine and hyperbolic functions in 
(O) can be expanded in the series. Hence, the function 
can be rewritten as 



* = 



«(k) = 



B z 



ik z V(k) 



l + \ 2 ab kl + \lkl 



(21) 



4vr C/ 44 + U m + (B2fc2/47T$ )*(k) ' 



(18) 



In order to describe the field distribution of a JV in 
the crossing lattice, the averaged magnetic induction Bj 
along the x-direction generated by both JV and inclined 
PV lines can be introduced. By substituting the found 
displacement ( [l8|) into equations (0,^), the magnetic in- 
duction of JV Bj = hj + h* is rewritten as 



Substituting this expression for 'J into (|T^,^) and omit- 
ting the difference between k and q, the equations for the 
dense PV lattice (which determine the field distribution 
and the energy of JV) are deduced as 



1/X 



' dq y 



l/Aj 



2tt 



l/s 



dq z 



l/s 



2tt 



dq z dq y <$> e lq y y+iq * z 
(2tt) 2 l + A ^ + A 2 b(? 2 



(22) 



Ei dkydk 



ik z u(k) 



and 



(2tt) 2 l + A 2 (fc,-Q,) 2 + fc 2 A 2 b 



i(k y 



„ 4> 2 f 1 ^' dq y f l ' s dq z 
oj = — 



(19) 



n t -l-i/Xj 2tt J_ 1/s 2tt 



where q y and q z are the wave vectors of a "bare" JV 
(\q y \ < 1/Aj, \q z \ < l/s), while the wave vectors k of the 
PV lattice are restricted also by the first Brillouin zone 
of the PV lattice (\k y \ < min(l/Aj, n/b), \k z \ < l/s). 
To get the energy £j it is necessary to add F cross to 
the energy of the JV itself. Obviously, the energy of a 
JV in the presence of PV lines is always lower than one 
of a "bare" JV. Indeed, for the displacement of PVs u 
determined by equation (p"8]), the energy F cross takes the 
minimum value which is smaller than zero, since F cross = 
at u = 0. Finally, the energy of JV in the crossing 
lattice obeys the equation 



1 + \lq 2 y + \\ b ql + B*q*/(MU 44 + C 66 q*)) ' 



(23) 



The last undefined parameter, Aj, can be obtained from 
the condition \dBj(y « Aj,z = 0)/dy\ ~ (47r/c)j c : 



IT fV X ' f 1 ^ 

-2- ~ Aj / dq y / dq z 

"V s J-l/A., J-l/s 



1 + XW y + \ 2 ab ql + qlBl/(An(U AA + C 66 g 2 )) ' 



(24) 



Bl 



dq y dq z 



1 



(2tt)2 1 + A^ + A^I 



32tt 2 



k^(k)^(-k) 



(2tt) 2 C/44 + U 66 + (A;2S 2 /47r$o)*(k) ■ 



(20) 



Equations (p^-p0[) together with the condition that the 
current density along the c-axis should be smaller than 
the maximum current density j c , determine completely 
the behavior of the PV lines and the JV. However, in 
further analysis it is convenient to investigate the dense 
(7s 3> a) and dilute (7s <C a) pancake lattices separately. 



II. JOSEPHSON VORTEX IN THE PRESENCE 
OF DENSE PV LATTICE 



The region of integration is shown in Fig. 2a. The 
rectangular domain of possible wave vectors replaces the 
usual elliptical one due to a peculiar core structure of 
JV. In anisotropic London model, the core of an ordi- 
nary vortex is defined by the elliptical stream liae of the 
persistent current having the depairing valuetS. How- 
ever, in our case the maximum value of q z is determined 
by the layerness of the medium while the largest value 
of q y is restricted by the Josephson critical current along 
the c-axis. The rectangular domain of wave vectors (Fig. 
2 a) can be divided into "screened" (\k z \ > 1/6) and 
"remote" (\k z \ < 1/6) subdomains. The first one corre- 
sponds to the region where one can roughly neglect the 
weak logarithmical dependence on k z in eq. (|l4j) to ex- 
press the tilt energy as: 



U44 ~ U44 



044k 2 , 



(25) 



For the case of the dense PV lattice, many PV rows 
are placed on the nonlinear JV core (Fig. lb, left sketch). 
It means that the magnetic field of a bare JV varies on 



with U44 = (B z $o/327r 2 A4 fc )ln(l + k z /(Xj + b~ 2 )), 



C 



- g»go In 
44 - 32^A 2 111 



f — =s 

lfc- 2 +(fc,-, 



/7 2 )+A<: 



and k z ~ yjl/bt 
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FIG. 2. The integration region in equations ( Jig| ) and ( |2c| ) 
for the dense PV lattice (a) and the dilute PV lattice (b). In 
the dense case (a), the region of the available wave vectors k 
of the PV lattice coincides with the accessible domain of wave 
vectors q; the continuous approximation (^) for <P is always 
valid. The approximation (E5I) is correct in the "screened" 
domain of wave vectors, while it fails in the "remote" re- 
gion. In the dilute case (b), the continuous approximation 
is correct only in the region \q v \ <C ir/b; \q z \ <C y/b; the 
"non-screened" region \q y \ > n/b is not accessible for k. The 
continuous approximation for Vl'(k) is broken in the "pinning" 
region \k y \ < ir/b; 7/6 < \k z \ < 1/s. 



16n 2 X e l f X 



' Cqq \ a b 

u^x 2 ab a7 



In- 



Xj 



(27) 



where 

A C ut ~ 



min(A c , min(76, max(6A c /A^{ / , \J bX c X ab / X e J b f ))) , 



while the renormalized penetration depth A a 
pressed as 



K ab 




B' 2 
A-kUaa 



(28) 



The physical reason of the renormalization of the in- plane 
penetration depth is related to the screening of the JV 
field by currents born by the local inclination of PV lines. 
From eq. ( p6| ) it is easy to see that the size of the non- 
linear JV core also decreases due to the interaction of 
JV and PVs, The similar conclusion was given earlier 
by KoshelevB, who considered the additional phase vari- 
ation of the order parameter born by the displacement 
of PVs. However, the shear contribution.. to the renor- 
malization of 8 j and Xj was neglected inQ, which could 
be done only for Xj > X a b (see equations ( pr^ ) and (p7|)). 
In the opposite case, i.e., when the London penetration 
depth exceeds the JV core size, the shear deformation 
becomes relevant and, as a result, Xj decreases with B z 
slower than it was proposed inQ. 

To understand how the field of JV is distributed in 
the real space, we rederive the results considering the 
free energy functional of the displacement u defined as 
a function of the spatial coordinates. In the limit 7s 3> 
a, the JV field varies on scales larger than the distance 
between PVs even near the JV core. This means that 
the field h p along the cc-axis can be averaged out on the 
scale larger than a: 



X , 



d 2 h; 



d 2 h* 
c dy 2 



B 



du 
dz 



(29) 



However, the short range variations of the field h p give 
the shear energy J7 66 and the tilt energy t/44. After ignor- 
ing the slow logarithmic dependence on k in the expres- 
sion for U44, one can conclude that the density of the tilt 
energy in the real space is U44U 2 (y, z), while the density 
of the shear energy is Uqq = CQ§{du/dy) 2 . Thus, the free 
energy functional is expressed as 



The expression (|25|) is completely wrong in the "remote" 
region in which it is necessary to use equation (|l5|). By 
using approximation ( Ba) for the tilt energy, the integrals 
in equations (|2^) and p4|) are easily evaluated (Appendix 
B) and we get 



X, 



X, 



max 



ye//' 



xeff 
^ab 



Xab 



(26) 



F cross = i J d 3 R^47rC 66 

+ AttUhu 2 + hlBz^T- + 2h ./ 

1 dz 



du 

dz 



(30) 



The first three terms represent the elastic energy (born 
by shear, electromagnetic tilt and Josephson coupling tilt 
rigidity, respectively), but the last term is related to the 
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interaction of the PV lines with the current generated by 
JV. 

In order to get the complete set of equations for the 
displacement u and the averaged magnetic induction Bj, 
we have minimized the functional (pQ) and have added 
together equations (|j) and (|2^): 



with the length A* b = . However, the set of equa- 



te. 



A 



d u 
dy 2 
d 2 B. } 
5 dz 2 



4ttU 44 u - 2B : 



dBj 

dz 



0, 



A: 



■ d 2 Bj 
dy 2 



B.£. (3D 
oz 



This set of equations is applicable if the continuous ap- 
proximation is valid (Aj ^> a) and, strictly speaking, only 
when the tilt energy Uaikz) can be replaced by the con- 
stant U44. The last condition fails on the distances far 
from JV (z 2 + y 2 /j 2 > b 2 ). In this "remote" region, the 
constant U44 has to be substituted by —C e JJd 2 jdz 2 . Be- 
sides, if X a b > 7S7 the parameter U44 should be replaced 
by — C44U 2 jdz 2 near the JV core (z < X a bh)- 

Even though we consider only the situation when the 
set of equation ( |3l| ) is valid, i.e., the case X a b < 7s and 
the region z 2 + yr/j 2 < b 2 , the solution of equations (|3l| ) 
seems to be quite complicated. The relation between the 
displacement u and the magnetic induction Bj, which is 
obtained from the first equation of |3l]), becomes nonlo- 
cal due to the shear rigidity of the PV lattice: 



B z 



8ir\/ CqqU 44 



dy 5 e 



\y-y\/s 



(32) 



where S — \ a b\J Cqq/{U 44\ 2 ab ) ~ A a 6 is the characteris- 
tic length of a nonlocality. However, the nonlocality is 
irrelevant if the space scale of the variation of Bj is sub- 
stantially large then 8, i.e., if Aj 3> A ab . In such a case, 
the equations ([n|) for Bj and u can be decoupled 

B z dBj 



AirU 44 dz 



Bj (Kl f ) 2 



d 2 Bj 
dz 2 



(33) 



The equation (|33| ) for induction Bj is the London equa- 
tion with the renormalized in-plane penetration depth 
X e JJ ' . Therefore, the field distribution Bj, not far from 
the center of the Josephson vortex (z 2 + y 2 /j 2 < b 2 ), can 
be approximated as 



B.j = 



-Kn 



\/z 2 /(\ e a f b f ) 2 + y 2 /\ 2 ) , (34) 



2-A^A_ 

where Ko(x) is a modified Bessel function of zero order. 
Using the free energy functional (|o[) and equations , 
it is easy to show that the energy of JV is determined by 
the field in its center, i.e., £j = &o/(8it)B j(y « Xj,z m 
8): 



£.1 = 



i^ 2 \T\ 



HKJs) 



(35) 



tions (31) becomes incorrect in the region z 2 + y 2 /j 2 > b 2 



which cuts off that length as A* b « b. Thus, the expres- 



sion (|3J]) coincides with the earlier obtained equation (|27 
in the studied case Xj > \ ab . The results ( |3^ , |35| ) can 
be interpreted in terms of the effective anisotropy pa- 
rameter 7 e -^ = X c /X e JJ which governs the JV lattice. 
Since X e J h > X a b, the effective anisotropy -f e ^ is re- 
duced in the presence of PVs with respect to the "bare" 
one 7 = XcJXab- The similar anisotropy 7 ej * was earlier 
introduced!] as a ratio — Xj/s, but these two dif- 
ferent definitions of 7^ give the same value in the case 
Aj > A a fc when the shear deformation is irrelevant. 

Here, we discuss how the core size and the JV en- 
ergy are changed with the magnetic induction B z if 
7s > A a ft. For quite high magnetic inductions B z > B\ = 
(&o/X 2 b ) x ('fs/Xab) 2 , the size of the nonlinear core Xj is 
smaller than X a b and the shear contribution to the free 
energy is important. The second logarithmic term in the 
JV energy (|27j) can be omitted, and the core size obeys 
the equation Xj(B z ) ~ y/jsa. With decreasing of induc- 
tion, the core size increases proportionally to B z and 
reaches X a b at B z « B\, At low fields, the shear interac- 
tion between rows is irrelevant, the JV core size becomes 
A,/ = X c s/X e JJ ~ jsa/Xab, and the energy of JV is deter- 
mined by the logarithmic term in (|27|), Below the field 
&o/X 2 b , at which the distance a between PVs exceeds 
X a b, the currents generated by PVs practically do not in- 
fluence on the JV field and, thus, the renormalization of 
X a b, A j and the £j vanishes. In the case of X a b > js 
the physical picture is different from the previous situa- 
tion. The core size obeys the law Xj ~ \Ffsa at fields 
B z > ®o/("fs) 2 - Below this field, the effective value of the 
in-plane London penetration depth X e JJ ~ X 2 b /a (jgjj ) is 
still larger than X a b, while the JV core size is saturated 
as Xj = js. This means that JV field shows different 
behavior far from JV (z 2 + y 2 /"f 2 > b 2 /j 2 ), where the 
redistribution due to the local inclination of PV lines is 
still important, and close to the JV core. 



III. JOSEPHSON VORTEX IN THE PRESENCE 
OF DILUTE PV LATTICE. 

Far from the JV center, z 2 + y 2 /^/ 2 > b 2 j^ 2 , the 
JV field varies slowly which causes the smooth varia- 
tion of the displacement u even for the case of the di- 
lute PV lattice (a > 7s). In that spatial region, the 
continuous approximation is still valid. On the other 
hand, near the JV core (\y\ < b), the JV current in- 
creases quite fast inducing a large displacement of the 
PV stack placed on the center of JV. In this case, the 
continuous approximation is not applicable. To de- 
scribe such physical situation, we consider the wave vec- 
tor area of k divided into two domains (Fig. 2b). In 
the first interval \k y \ < ir/b and \k z \ < 7/6, the func- 
tion \P can be still roughly approximated by equation 
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* « $ /(l + A?fc2 + A 2 b fc 2 ), while * « $ b/(2\ c \ ab k z ) 
in the second region \k y \ < n/b and 7/6 < \k z \ < 1/s 
("pinning" region in Fig. 2b). Following this approach, 
the energy of JV is evaluated as 

1 

X 1 + \ 2 c q 2 v + \\ h ql + B*q*/(**(UM + U m )) 

, *8 J* 



167r 2 A c A afc \7s 



128^ 3 aA 2 A 2 b 7 7/b U 44 + B z <S> k z /{8nX c X ab a) ' 



(36) 



The first term comes from the spatial region far from the 
center of JV while the second and the third terms are 
related to the vicinity of the JV center. The screening 
of the "bare" J V field vanishes near J V ( "non-screened" 
region in k-space) which determines the second term in 
(|36|). The last term in ( |36| ) represents the energy gain 
due to the strong interaction between the PV line placed 
on the JV core and the JV currents (the energy gain 
of a PV stack placed on a JV in the limit 7s ^> X a b 
and B z — > was calculated by KoshelevQ). Since the 
last term is sensitive to the mutual position of JV and 
the nearest PV line, this contribution can be called as 
the "crossing lattice pinning" . Using the results of Ap- 
pendix B and taking into account that the evaluation 
B z &ok z /(8A c A a ba) ^ Cakl is held in the "pinning re- 
gion" (k z > 7/6), the energy of JV is finally obtained: 



£j 



$1 ( 2^iC 44 7 2 / t4 c M Kb 



$5 



X C ut 



In 



^ m ( A 

167r 2 A c A a & \7s 



47raA, 



■ arctan 



b — 7s 



'Un/C^sh + 7W C44/U44 



, (37) 



where fi = B z <I>o/(327r 2 A c A 2 b \/ C44U44) < 1 is the di- 
mensionless function depending quite slowly on B z and 
the numerical parameters \i\ and fi2 are about unity. 

Next, we will discuss how the renormalization of the 
JV energy comes in with increasing of the z-component 
of the magnetic field. At low fields, B z <c ^o/^ab 
(a 3> X a b), the first term and the last term in d37| ) can 
be omitted and the expression for the energy of a "bare" 
JV reported earlier int3c3 is reproduced 



S.T = 



16n 2 X c X a b 



In |'*£ 

7s 



(38) 



For the case X a b > is, the renormalization of JV en- 
ergy becomes relevant at B z w <I>o/A 2 b , i.e., earlier than 
the JV core size starts to decrease which occurs only in 



fields B z > <E>o/(7s) 2 . The origin of this behavior is that 
the additional current along the c-axis induced by tilted 
PV stacks is much smaller than j c near JV core in the 
field interval $o/Xl b < B z < ^ /(js) 2 , but the inclina- 
tion of all PV lines can still cause the renormalization 
of the JV field on scales larger than a. In the field in- 
terval <$>o/X 2 ab < B z < $ /7 2 s 2 ^the main contribution 
to the Josephson vortex energy ( p7|) comes from the first 
term related to the tilt elastic rigidity (born by Josephson 
coupling of PVs) and shear elasticity of the PV lattice. 
Strictly speaking, from our rough estimation of (pp[), we 
can not conclude how strongly £j is suppressed in that 
field interval, i.e. in the presence of the dilute PV lat- 
tice. Nevertheless, the pinning energy (last term in fl37|)) 
could be the same order of magnitude as the first and 
the second terms in ( |37| ) in fields B z ~ ^o/^ab an d ma y 
decrease £j substantially. 

Another interesting possibility arising due to the 
"crossing lattice pinning" is the rearrangement of the 
PV lattice in the presence of the JV sublattice. In the 
in-plane magnetic fields B x , JVs form a triangular lat- 
tice with distances aj and bj between JVs (see inset in 
Fig. 3a). In general, the PV sublattice and the JV sub- 
lattice are not commensurate: a,j pb with integer p. 
This means that the considered one-component displace- 
ment of PVs u = (u(y, z),0, 0) (the "shifted" PV lattice 
shown in Fig. 3a) does not provide the energy gain com- 
ing from the "crossing lattice pinning" since the PV rows 
can not occupy the centers of JVs. However, the PVs can 
be rearranged in order to occupy all JVs (the "trapped" 
PV lattice shown in Fig. 3b) if the PV lines shift also 
along the y-direction: u = (u x (x, y, z), u y (x, y, z), 0). 
The "crossing lattice pinning" decreases the free energy 
of the "trapped" PV lattice, while the additional shear 
deformation acts in the opposite way through increasing 
the free energy. For the case B x < jB z , the energy gain 
related to the "trapped" PV lattice is calculated by nor- 
malizing the last term of equation per unit volume: 



B x $ 

E tr = \x — arctan 

AnaXr 



b — 7s 



'V 44/C '44S& + J\/ C44/U44 



(39) 



But, in order to trap the PV lattice, the total displace- 
ment of PVs along the y-axis between the two nearest JV 
rows, i.e., on the scale aj, should be about b. Following 
the simple analysis 2 ^, the extra shear deformation (inset 
in Fig. 3b) is about Sb/b ~ b/aj (6b is the change of 
the distance between rows of PVs) and the energy loss 



can be estimated as: 



E 



shear 



aj 



B x 

Wz 



(40) 



with numerical constant v < 1. For the case 7s X a b, 
the shear elastic energy ( ^0| ) is strongly suppressed in the 
fields B z < $o/(7s) 2 where the "crossing lattice pinning" 
is active, since Cqq is exponentially small if a > X a b (13). 
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FIG. 3. The different substructures of the crossing lattice: 
a) the "shifted" PV lattice characterized by one-component 
displacement along JVs, b) the PV lattice trapped by the 
JV sublattice for the case when the distance aj between JVs 
exceeds the distance b between PV rows, c) the "trapped" 
PV lattice for a.j < b, i.e., in the case of the field orienta- 
tions very close to the afo-plane. The dotted lines depict the 
rows of unperturbed lattice (crossings of these lines and filled 
circles mark the positions of unshifted PVs) in all sketches. 
Dashed-dotted lines indicate the rows of the "trapped" PV 
lattice which are deformed in order to match with JV sublat- 
tice. The arrows directed from the filled circles to the open 
ones show the two component displacement of PVs for the 
cases represented in sketches b) and c). Inset in a): the JV 
sublattice with lattice parameters aj and bj (lines mark the 
Cu02 planes). Inset in b): the additional deformation of the 
PV lattice which is required for trapping of PVs by JVs. The 
upper sketch is the equilateral triangle of the unperturbed PV 
lattice which is incommensurate with JV lattice (aj 7^ pb with 
integer p) . The lower sketch is the isosceles triangle of the P V 
lattice matched with the JV sublattice (oj = p(b + 5b)). 



Therefore, the "trapped" PV lattice seems to be realized 
as soon as a >r-^s. In the opposite case, A h > 7s, the 
transformationE2l from the "shifted" PV lattice to the 
"trapped" PV lattice occurs when the energy gain E tr 
exceeds the energy loss E s h ear . It happens in a certain 
out-of-plane field between the field <&o/(7s) 2 , at which 
the "crossing lattice pinning" is activated, and the field 
B z ~ &o/^ab' where the shear elastic energy rapidly de- 
creases. Next, we discuss the difference between the con- 
sidered "trapped" state and the "chain" state proposed 
for the crossing latticed. The "trapped" state is related 
to the rearrangement of PVs on the scale aj between the 
nearest rows of JVs. On the other hand, the "chain" 
state is associated with the creation of an extra PV row 
(an interstitial in the PV lattice) on a JV, but the influ- 
ence of the neighbouring JVs is completely ignored. As 
a result, the "trapped" and "chain" states have the dif- 
ferent in-plane field dependence of the out-of-plane tran- 
sition fields. The out-of-plane transition ficldEil between 
the "shifted" and "trapped" PV lattices does not depend 
on H ao in contrast to the i? a (,-dependent out-of-plane 
fieldLJ of-,the destruction of the "chain" state. Since the 
analysisQ is correct only in the case of 7s 3> A a b and 
o. ^> A a b, the transformation of the PV lattice discussed 
here seems to be more likely in the case X a b > 7s. 



IV. PHASE DIAGRAM OF VORTEX LATTICE IN 
TILTED MAGNETIC FIELDS 

In this section we discuss the vortex lattice structures 
formed at different field orientations. The tilted lattice 
consists of mono-oriented vortices and transforms contin- 
uously from the tilted PV stacks in fields near the c-axis 
(Fig. 4a) to the long JV strings connected by PV kinks 
for the field orientations close to the a6-plane (Fig. 4b). 
On the other hand, the tilted lattice is topologically dif- 
ferent from the crossing vortex structure (Fig.. 4c), and 
they replace each other via phase transition!!. For the 
analysis of the vortex phase diagram in tilted fields, the 
free energy of the crossing and tilted vortex structures 
will be compared. We concentrate .an the case 7s X a t, 
when, according to Bulaevskii et aln and KoshelevEl, the 
tilted lattice-.is energetically preferable above the lock- 
in transitior£l. We will consider a thin superconducting 
platelet with the c-axis perpendicular to the plate. In 
this geometry the lock-in transition occurs at very low 
fields! B,»(l- 7j 2 )$o/(47rA2 b )ln(7s/£ Qb ) with demag- 
netization factor n z (l-n z <l). 

For the field oriented close enough to the c-axis, tan 9 = 
B x /B z <C 7, the free energy of the tilted lattice F t can be 
evaluated as F t = if + \C^{^ = 0)Bl/B 2 z in analogy 
to the analysis given in ref.Q. Here, represents the 
free energy in the absence of the in-plane magnetic field, 
while the tilt modulus is expressed as (^"(k = 0) = 
B 2 JA-k + C e A { S with Cl{ f defined in (H) for the case of 
B z > ^o/(47rA^ fc ). As a result, we have: 
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tice while the interaction of PVs and JVs is taken into 
account through the renormalization of the JV energy: 




FIG. 4. The 3D sketches of the different vortex structures 
in the tilted magnetic field with the components H c and H a t 
along the c-axis and in the afe-plane, respectively: a) the tilted 
vortex lattice near the c-axis (TI), when the current between 
CUO2 planes is much smaller than the critical value j c , i. e., 
the Josephson strings linking PVs are not developed; b) the 
tilted vortex lattice far away from the c-axis (Til), when the 
JV strings are formed; c) the crossing vortex lattice. 



F t ~ -2- 



8tt 



+ 3.68 



H, 



B 



8tt 



Bl 



o 



hi- 



H, 



c2± 



2(47rA afc ) 4 Bl 64ir 2 \ 2 B z B z 



(41) 



where H c2 ± = $o/2vr^ fc . The first two terms form the 
free energy for B x = 0. The third term is the in-plane 
magnetic energy, the fourth one comes from the electro- 
magnetic interaction of the inclined PVs, and the last 
contribution is connected with the Josephson coupling of 
PVs. 

The free energy of the crossing lattice F c consists of two 
contributions from the PV sublattice and the JV sublat- 



F ^B1 

c ~ 32tt 2 A 2 , 



<f> B z 



In- 



B, 



^ + ^£j. 

8?r $ 



(42) 



The renormalized JV energy, Ej, is defined by equa- 
tion (^o|) in which the lower limits of integrations are 
restricted by the conditions q y , k y > 1/aj and q z , k z > 
l/bj. " 

The tilted lattice is energetically preferable in the fields 
oriented near the c-axis because F t oc _B 2 , while F c oc B x , 
s., Ft < F c for low B x . The phase boundary between 
the tilted lattice and the crossing structure can be ob- 
tained from the condition Ft — F c which is rewritten in 
the form: 



B,, 



£j 



Bl 



$ 1.84$2/( 47 rA ab ) 4 



$ B z /{6An 2 \ 2 c )HH c2± /B z y 
(43) 



The transition from the tilted lattice to the crossing 
structure occurs at the field oriented quite close to the 
c-axis for high anisotropic superconductors due to: 1) 
the high energy cost of the inclination of PV stacks in 
the tilted lattice related to the electromagnetic interac- 
tion of PVs, and 2) the decrease of the JV energy in 
the crossing lattice structure. For the dense PV lattice 
B z 3> <&o/(7s) 2 and A a fc > 7s, equation ( [l3| ) can be sim- 
plified: 



B, 



a 



66 



B 2 



7$o 



4.3tt 2 \ 2 



^\n(H c2 jB z ) 



(44) 



Next, we will study the field orientations close to the 
a6-plane, B x > "fB z . Here, the electromagnetic interac- 
tion between PVs in the tilted lattice is not so important 
and the fr/ae energy in the low c-axis fields B z < $0/^ 
is reducedQ to: 



$oB x 



■In- 



87T 327r 2 A a6 A c -ys 2 B x 



H,jB z 
An 



(45) 



where Hj = $o/(47rA 2 h ) ln(7s/£ a (,). The first two terms 
are related to the energy of JV strings while the last one 
is associated with the energy cost of the formation of 
PV kinks. In more detail, the PV kink generates the in- 
plane current which decreases with the distance r from 
the kink center as 1/r up to a critical radius r$ of the 
region with 2D behavior where the current along the 
c-axis is about the maximum possible current j c u. At 
larger distances, the in-plane current decays exponen- 
tially. Simple evaluation gives tq = 7s for the PV kinkQ. 
We note, that the tilted vortex lattice in the considered 
angular range of the magnetic field orientations seems to 
exist as a kink-walls substructure, where kinks (belong- 
ing to different vortices) are collected in separated walls 
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parallel to the yz-planeB For the kink-wall substruc- 
ture of the tilted lattice, the contribution to the free en- 
ergy Jl5|), attributed to the PV kinks, is slightly reduced 
in the high in-plane magnetic fields B x > <&o/js 2 ni, 
which can be taken into account through renormaliza- 
tion Hj = $o/(87rA2 6 )ln( 7 if c2 /S x ). 

In the considered field interval, B x > jB z , B z <C 
$o/\ 2 b , the renormalization of the JV energy in the 
crossing lattice structure vanishes. However, the inter- 
action of PV and JV sublattices still manifests itself 
through the "crossing lattice pinning" : 



F = ^ 
c 8ir ' 32^ 2 A ab A c "* ~/s 2 B x 



■ In ■ 



H c i±B z 

ill 



[iB z 



1 gxfo 



■ arctan 



1 — \J B x /H a 



, (46) 



with H clJ _ = $ /(47rA 2 h )ln(A afc /^ Q t,) (the critical radius 
ro for the in-plane currents of a PV stack is about X a m) , 
H a = $o/js 2 and H% = 7<& /A 2 b . The third term is the 
energy of the unperturbed PV lattice, while the last term 
corresponds to the "crossing lattice pinning" contribution 
which can significantly decrease the free energy F c in the 
in-plane field interval 7$oAa& < B x < $ /js 2 . The 
difference between the "crossing lattice pinning" contri- 
butions to the free energy in the cases of low B x < yB z 
and high B x > "fB z in-plane fields (see equations pS|), 
(|46|)), emerges because the number of PV lines is suffi- 
cient to occupy all JVs (Fig. 3b) at B x < jB z while 
some JV strings do not carry PV rows (Fig. 3c) in the 
opposite case. By analyzing equations ([15]) and fl46|), 
we can conclude that, at least for B x < 7 < &oMofc anc ^ 
B x > $o/7s 2 , the tilted lattice exists near the a6-plane 
since the condition F t < F c is held due to the inequal- 
ity Hj < H c i±. The tilted lattice is replaced by the 
crossing lattice with increasing the out-of-plane magnetic 
field above B z — B x /j. However, it is difficult to deter- 
mine the contour of the possible phase line between the 
crossing and tilted vortex structures, since it requires the 
more precise calculations of the free energies F c and F t 
in the region B x < 7-B z . In the intermediate in- plane 
magnetic fields 7^o/A^ 6 < B x < <fr /-fs 2 , the "crossing 
lattice pinning" could make the crossing structure to be 
more energetically preferable with respect to the tilted 
lattice. In that case, the crossing lattice (CII, Fig. 3b) 
with a < aj transforms into the crossing lattice (CIII, 
Fig. 3c) with the extremely dilute PV sublattice a > a.j 
at the angle 9 = arctan B x /B z ~ arctan 7. 

Therefore, we find a complicated picture of phase tran- 
sitions between the tilted vortex structure and the cross- 
ing vortex structure in the case 7s < A Q (,. The proposed 
phase diagramed is shown in figure 5. As it was suggested 
earlieiu and according to our calculations by using equa- 
tion (f43|), the tilted vortex structure (TI) of inclined PV 
stacks (see Fig. 4a) can be replaced by the crossing lat- 
tice quite close to the c-axis (see phase diagram obtained 
for 7 = 500, inset in Fig. 5). 
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FIG. 5. The proposed phase diagram of the vortex solid 
phase in the oblique magnetic fields calculated using equations 
©. ©1 @j @ and © with parameters mentioned in 
the text, T = 45 K, 7 = 100 and v = 1. The dotted line is 
the line B x — jB z , while the shaded area marks the region 
inside which the transition from the crossing lattice C(II) to 
the tilted lattice T(II) or the crossing lattice C(III) happens. 
The arrows from enframed TI and Til are directed toward the 
regions where these vortex structures are realized. Inset: the 
part of the phase diagram close to the c-axis for strongly ani- 
sotropic superconductors with 7 = 500 (7s > \ a b) and T = 45 
K; the solid line marking the transition from the tilted lattice 
(TI) to the crossing lattice is obtained from eq. (^), while 
the dashed line, corresponding *o the same transition, is cal- 
culated by using eq. (6) of Ref.u. The parameters are chosen 
to give some insight to the behavior of Bi2Sr2CaCu20s-i in 
the oblique magnetic fields. 



Nevertheless, in the high c-axis magnetic fields, the calcu- 
lated in-plane magnetic fields of this transition are higher 
than ones obtained by using the modelQ. This difference 
comes from the shear contribution to £j which was omit- 
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ted in reffl. The crossing lattice, which exists in a wide 
angular range, can have different substructures. At high 
enough out-of-plane fields B z > $ /7 2 s 2 , the "shifted" 
PV sublattice is realized in the crossing lattice structure 
(CI, Fig. 3a). In this substructure, the JV currents shift 
the PVs mostly along the x-axis. The "shifted" phase can 
transform into the "trapped" PV lattice (CII, Fig. 3b), 
when the energy gain related to the "crossing lattice pin- 
ning" exceeds the energy needed for the additional shear 
deformation (the dashed line in Fig. 5 separating CI and 
CII has been obtained from the condition E tr — E s h ear ). 
Around the line B x — ^B z , the lattice CII can be changed 
by the tilted lattice (Til) with JV strings linked by PV 
kinks (Fig. 4b) or by the crossing lattice structure (CHI, 
Fig. 3c) at which all PV stacks are placed on a few JVs. 
The domain in the H c — H a b phase diagram with the lat- 
tice CIII is determined by the condition F c < F t where 
F t and F c are defined by the equations ( fi"5| ) and (|4^), 
respectively. With increasing temperature, the region of 
the lattice CIII becomes narrower and disappears at a 
certain temperature (see Fig. 6). The proposed phase 
diagram suggests the possibility of the re-entrant tilted- 
crossing-tilted phase transition as the magnetic field (at 
least with low B < 7$ M 2 h or high B > $o/7s 2 abso- 
lute value) is tilted away from the c-axis to the a6-plane. 
Such possibility for low fields was earlier mentioned in 
the worksU'ES in which the interaction of crossed sub- 
lattices was not considered. Moreover, we note that the 
instability of the tilted-iattice was found numerically by 
Thompson and MooreO (for 7 < 100) only at the inter- 
mediate field orientations 0° < 6\ < 9 < 8 2 < 90°, which 
could also support the discussed scenario. The parame- 
ters taken for the H c — H a i, phase diagram at T = 45 K 
(Fig. 5) and the H a b — T phase diagram at the magnetic 
field orientation B z = B x /j (Fig. 6) were chosen to give 
some insight into the behavior of vortex array in BSCCO 
in the tilted magnetic fields (see further discussion) as 
A ab - 2000/0 -T 2 /T c 2 A, U = 30/0 - T/T c A, 
s = 15 A, T c = 90 K, 7 = 100 (7 = 500 for inset in 
Fig. 5 and 7 = 150 for inset in Fig. 6). 



pearance of the hexagonal-order along the c-axis found 
by neutron measurements^. With further tilting of the 
magnetic field, the linear dependence of the c-axis melt- 
ing field component Ji" 1 ^^) abruptly transforms into a 
weak clepcndenccEZrEj, which, as was shownES, can not be 
explained in the frame of the modeD. Such behavior sug- 
gests a phase transition in the vortex solid phase in tilted 
magnetic fields in Bi2Sr2CaCu20s+<5, which w.as .detected 
in the recent ac magnetization measurementscJ'Eil. 
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V. CONCLUSION 



This theoretical investigation was partially motivated 
by the recent intensive experimental studies of the vor- 
tex lattice melting transitipJpEj as well as transitions 
in the vortex solid phascc3 : Eil in Bi2Sr2CaCu208+a sin- 
gle crystals. The observed linear decay of t he .. Cj axis 
melting field component if™ with in-plane fieldsESEZl was 
interpreted!] as an indication of the crossing vortex lat- 
tice. Thus, the tilted lattice could be replaced by the 
crossing vortex structure quite near the c-axis. Accord- 
ing to our calculations, the angle where such transition 
may occur is about 7° at B z ss 100 Oe and 7 » 500 
while that angle reaches 14.5° in the higher out-of-plane 
field B z — 500 Oe, which correlates well with the disap- 



FIG. 6. The H a b — T phase diagram of the vor- 
tex solid phase at the field oriented near the afe-plane 
(B x /B z = 7 = 100). The region of the crossing lattice 
phase (CIII) becomes narrower and finally disappears with 
increasing temperature. Inset: the same phase diagram for 
7 = f50. The dotted lines correspond to the experimen- 
tally found temperature dependence H^l cx 1 — T 2 /T% of the 
in-plane characteristic fields of the vortex lattiee melting tran- 
sition in Bi2Sr2CaCu2 0s+a in the tilted fieldsEj (for instance, 
symbols represented in the inset exhibit the temperature de- 
pendence of the maximum in-plane field of the vortex lattice 
melting transition). 
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As was mentioned by Ooi et a/0, the behavior of the new 
anomaly of the local magnetization in BSCCO attributed 
to the phase transition in the vortex solid slightly reminds 
the peak effect related to the vortex pinning, which, in 
turn, could be induced by the vortex trapping by planar 
defectsEJ. Such analogy, as well as a very weak depen- 
dence of the c-axis magnetic field component H* Qf-jthe 
novel phase transition on in-plane magnetic fieldsEU in 
a wide angular range, could suggest the tansition from 
the "shifted" PV sublattice (CI) to the "trapped" PV 
sublattice (CII) in the crossing lattice structure. Our es- 
timation of the c-axis field H* rap of the transition from 
CI to CH, (H^p fa 450 Oe for 7 = 100 at T = 45 K (see 
Fig. 5)) — is in a reasonable agreement with experimental 
findings^ H*(T = 45 K) ps 430 Oe. Near the a6-plane, 
the properties of the observed anomaly-changes abruptly 
and the field H* sharply-, goes to zeroEil, which may in- 
dicate a transformations in the JV lattice or a trace of 
the phase transition from Til to CILL At higher in-plane 
fields, the found step-wise behaviorEa of the vortex lat- 
tice melting transition may be related to the existence of 
one more phase transition in the solid phase. Interest- 
ingly enough, all characteristic in-plane fields of the vcp=j 
tex lattice melting transition depend on temperatureEJ 
proportionally to 1 — T 2 /T 2 , which is similar to the cal- 
culated temperature dependence of the phase transition 
from CIII to Til (see Fig. 6). 

In summary, we discussed the crossing lattice struc- 
ture in a strongly anisotropic layered superconductor in 
the framework of the extended anisotropic London the- 
ory. The renormalization of the JV energy in the crossing 
lattice structure was calculated in the cases of the dense 
PV lattice as well as the dilute PV lattice. It was shown, 
that the "crossing lattice pinning" can induce the rear- 
rangement of the PV sublattice in the crossing lattice 
structure as soon as the out-of-plane magnetic field be- 
comes lower than a certain critical value. The free energy 
analysis indicates a possibility of the re-entrant tilted- 
crossing-tilted lattice phase transition with inclination of 
the magnetic field away from the c-axis to the a6-plane 
in the case of X a b > 7s. 



with real numbers r, t, and a. 

Next the sum ^i{k y ,k z ) = *oEQ B=2ra /(,/( fc z>^ - 

Qy)/0- + ^ab^z + ^c(K ~ Qy) 2 ) needs to be estimated. 
By using inequality k z < 1/s, we obtain f(k z , k y — Q y ) ps 
fj(k y — Q y ) where fj ps 1 for \k y — Q y \ < 1/A,/ and zero 
otherwise. In the case of the dense PV lattice (6 <C Xj) 
we retain only the term with n = in the sum and get 
the expression ^i(k y > 1/Aj) = and ^i(k y < 1/Aj) = 
^o/(l+X 2 b k 2 +Xlk 2 ). Taking into account the inequality 

k y < 1/Aj < l/b and yjl + X 2 ab k 2 b/X c < 6/(73) < 1, 
one can rewrite ^i(k y < l/Xj,b -C A,/) ~ ^{k y ,k z ). 
Thus, we come to the equation (|16|). 

For the case of the dilute PV lattice (Aj <C b), 
many terms give contributions to the sum ^ ~ 
*o J2n=~N 1 /( 1 +Aa6 fc ?+Ac(fcy-27rn/6) 2 ), since inequal- 
ity \k y — 2-Kn/b\ < 1/Xj is held until n exceeds N ^> 1. 
Thus, the function \I>i is estimated as: 



dx 



7T Jl/Xj 1 + A a 



L.2 

b^z 



X 2 x 2 



(A2) 



Finally, we obtain *r = (1 - (3(k y ,k z ))^> with 



P(kyi k z ) 



sMy/l + A^fc^/A^-cosfeb) 
sinh(yTTApf6/A c ) 



1 arctan 

7T 



A.7yr+Apf, 



(A3) 



At \k z \ -C 1/s, it is easy to show that (3 < Xj/b -C 1 and 
the approximation \&i = W used in ( |l6| ) is excellent. Only 
for k z ~ 1/s, the function ^ can differ from the function 
^ by a factor about unity in the case of the dilute PV 
lattice (the factor is 1 — f3 ps 0.5 in the framework of our 
rough consideration) . However, the correct estimation of 
the value of the factor depends on the type of the smooth- 
ing function and requires more precise analisys than one 
in the framework of the London approach. Therefore, 
we can always assume ^1 = W in our semiquantitative 
consideration. 



APPENDIX A. APPROXIMATE SUMMATION IN 
EQUATION (0) 

Our aim is to sum the sequence: 
^ J{k z ,k y - 2rrn/b) - (B z /2<P )ik z u(-k) 

*° E lkMk) — i + A^ + m-2w>) 2 — 



n— — 00 



APPENDIX B. EVALUATION OF INTEGRALS 
IN EQUATIONS (EHE3J36I) 



ik z ^x(k y , k z )u(k) + (B,/2$„)*2*(*v> *«)«(k)«(-k). (Al) 

The equation ( pj| ) for ^(k y , k z ) can be directly obtained 
by using the well-known mathematical equality 



E 



a sinh(ra) 



In this appendix, the integrals in equations ( |23| , |24 36) 
are evaluated. We start with the dense PV lattice, a <C 
Xj. In the region q z < l/b, the tilt energy is sma ll (|15|) . 
Therefore, the denominators of the integrands in (|23j2J) 
are substantially larger than ones in the case q z > l/b 
and we can roughly neglect the contribution related to 
the region q z < l/b. Using equation (p5) for the tilt 
energy in the domain q z > l/b, the integra s (|23|,|24|) can 
be rewritten as follows: 



r 2 + (t + 2im/a) 2 2r cosh(ra) — cos(ta) 
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1/Xj , r 1/s , m) 



dq y / dq 



l/Xj Ji/b l + X 2 c ql + X 2 ab q 2 + B 2 q 2 /(4ir(U 44 + C 4A q 2 + C 66 q 2 )) 

(Bl) 

where /(??) = <&q/327t 3 for (p3|) and / = Xjq 2 for (p3). After multiplying numerator and denominator by factor of 



E/44 + 6*44^ + C66?y, the expression (|Bl| ) is reduced to 

■ * * r m, jg«±g^t^MM) = (B2) 

The term (1 + X 2 q 2 + X 2 ah q 2 z )C 44 can be neglected with respect to B^/Att in the denominator. Indeed, the maximum 
value of X^qyCa (see eq. (p5|)) is about X 2 C 44 \/X 2 ~ B z $q/X % j <C Bf, while the maximum value of X 2 ah q 2 z C 44 is 
A^C44l/Aj x (Xj/js) 2 which is even smaller, since Xj < 7s (see eq. @)). Next, the integration in (B2) over q z can 
be taken easily: 



iA« V B '/ 47r + A '^ C/44 + C66(? y) s Ac(fe v /s2/4^ + A2 fc (Z744 + C 66g 2)/ 

where the function 0*, defined as 

0*(q v ) = arctan ( ^jsjl + £f /^ttA^m + C^gg)) J _ arctan ( ^- ^/l + Bf/ (47^(^44 + 0,6^)) J , 



determines the lower cutting value l/A CM t of q y . Namely, by using (|26|) we can assume that the argument of the first 
arc tangent in the last expression is larger than 1, for most of values of q y . Then, the value of 0*(q y ) is about 7r/2 if 
the argument of the second arc tangent is smaller than 1; otherwise 0* is close to zero. Thus, we obtain the expression 
for X cut - 



1/2 



min(A c ,min( 7 6,max(feA c /A:{^ V / &A c <5/A:{^))), 



where we denote 5 = W Cqq/U 44 ~ A a fc and take into account that X cu t should be smaller than A c . As a result, the 
integral QB3|) can be evaluated as: 



I = 4 / dq y f(q 2 y ) 4 4 - + 2^ / dqyfiql)— V (B4) 



1/Ae " B2/47T + A2 b ([/ 4 4 + CW ? 2).S 7 1/Aeu4 » A c g y ^ B 2/ 47r + + ^2) 

In case of the dense PV lattice, the first integral is small and can be neglected. In addition, we also can omit the term 
C 66 q 2 in the expression B 2 /4n + X 2 b C 66 q 2 . Finally, we have 



/(a<Aj) = 27T - / v - gg. B5 
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After taking /(<j^) = Xjq 2 in JB5| ) and by ignor- 
ing Ceeqy/U 44 in the case Xj > \ a b or by neglect- 
ing unity in the opposite case, we obtain the expres- 

A s 



sion Aj » max ^^J^ttj v-^b / ' ^ coinc i°- es w ith 

(^6| ) since (5 ?s A a b. The energy of JV (|^) is eas- 
ily derived if one puts f{ql) = $§/(32tt 3 ) in ®. 



Next, we roughly estimate the first integral in equation 
@. The inequality (1 + X 2 c q 2 y + X 2 ab q 2 z )C i4 < #J is 
still correct in the domain q z <C 7/6, q y <C 1/b for 
£?z ^> ^o/Ag. Thus, we can get the estimation (B4) 
also for the dilute case (a > 7 s). However, in contrary 
to the dense PV lattice, the contribution related to the 



first term in (B4) remains important: 



I(a > Aj) = 4= 



C447 



2tt 



1/Ac 



dq v f{q 2 v ) 



A \ eff 

A c A ab 



dq y 




S^flg- ( B6 ) 

(744 



Here, we have introduced numerical parameters and 
/U2, since the upper limits of integration are not well de- 
fined. The corresponding contribution to the energy is 
obtained from the last equation for / = $q/327t 3 as pre- 
sented in the text. 
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